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n^ ', A method is presented which ahows for a tremendous speed-up of computer 

C I simulations of statistical systems by orders of magnitude. This speed-up is 

Q ' achieved by means of a new observable, while the algorithm of the simulation 

remains unchanged. 
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I. INTRODUCTION 

In performing a Monte Carlo simulation, one or several observables are chosen, for which 
a simulation average is recorded. A common choice for such an observable is related to a 
histogram as introduced by Salsburg |jl]] and popularized by Swendsen and Ferrenberg P|, 
0]. The histogram allows for the estimation of the density of states, from which a variety 
of physically interesting system properties can be computed. Recently, Oliveira et al. @] 
suggested a new method, called transition observable method throughout this paper, which 
likewise enables an estimation of the density of states but leads to considerable reduction of 
the computing time. In other words, the simulation may be dramatically accelerated without 
modifying the algorithm but by changing the observable recorded during the simulation.0 

If the standard histogram method is applied to magnetic systems, the number of mi- 
crostates with energy E and magnetization M is counted during the course of simulation, 
i.e., every microstate yields one entry in an energy- magnetization histogram. From this 
histogram the density of states can be computed (cf. Sec. [11 C|) . 



In the transition observable method, each microstate of the Monte Carlo sample is exploited 
in a much more sophisticated way. For the particular realization of the method introduced 



in Sec. pT| , this means: In an extended histogram, the number of possible transitions is 
recorded from particular microstates (in the Monte Carlo sample) with energy E and mag- 
netization M to "neighbouring" microstates (not necessarily in the Monte Carlo sample) 
with energy E ± AE and magnetization M ± AM, which can be reached from the particular 
microstates of the sample by applying single spin flip operations. Again, from this extended 



histogram, the density of states can be computed (cf. Sec. |IIDD . 

As already pointed out by Oliveira , the advantage of the transition observable method 
is that every microstate of the Monte Carlo sample is investigated much more extensively 
than by the standard histogram technique. Therefore, given a certain sample of microstates, 
the density of states can be calculated more accurately from simulation averages of the 
transition observable introduced below than by standard methods. Additionally, as it is 
the selection of microstates using pseudo random numbers which is costly in computing 
time, the increase in computing time from such a more extensive exploration of the chosen 
microstates is absolutely neghgibleQ. However, the effect on the data quality is significant 
and can amount to orders of magnitude. In the case of the examples studied in this paper, 
we find an efficiency gain of roughly two orders of magnitude! This efficiency gain can be 
expected to grow proportional to L*^/^, the square root of the volume of the system. 

It is this enormous gain of efficiency which should motivate the reader to occupy himself 
with the underlying formalism, which is indeed simple to implement in a simulation, but is 
somewhat heavy to formalize. 



Section |IIA| and [II B| aim to familiarize the reader with the language used throughout this 



paper and with some aspects of the Monte Carlo procedure. In section [11 Cj, the standard 



^ Since this new observable represents properties of the system under consideration, it is of course 
completely independent of the particular simulation setup. 

^At least in the case of the particular realization of the method introduced below. 



histogram technique is reviewed in the context of the calculation of the density of states. The 



transition observable is introduced in Sec. IIP. In Sec. HE, it is shown that this observable 



includes the one presented in reference as a special case. The rest of the paper (section 



np is devoted to the comparison of the efficiency of computer simulations using the new 
method in contrast to a standard histogram technique. This is done for the examples of 2d- 
and 3(i-Ising systems where we find a speed-up factor of ~ 40 in the 32^ Ising system and 
^ 250 in the 10^ Ising system. 

II. CALCULATING THE DENSITY OF STATES BY MONTE-CARLO 

SIMULATION 

A. Conventions and notation 

In this paper, we use the language of discrete Ising systems on hypercubic lattices of 
linear size L m d spatial dimensions with Hamiltonian 

K(S') := —J 2^ ^i^i " h 2^ CTj 

= : E{S) - hM{S) , ^ e r^d , (1) 

where h denotes an external magnetic field. E{S) is the interaction energy and M{S) the 
magnetization of the particular microstate S* = cti, o"2, ..., ctj, ..., a^a (= particular configura- 
tion of the spins cxj, i = 1, 2, .., L"' on the L'^ lattice) with ai G { — 1,+!}. The configuration 
space of the Ising system is denoted by F^d, and {i,j) indicates a summation over all pairs 
of nearest neighbours. 

The discreteness of the Ising systems gives rise to a minimal energy and magnetization 
spacing, denoted by AE and AM, respectively. Summations over interaction energy E 
(magnetization M) cover all energy (magnetization) values accessible. 
In general, in order to simplify the notation, system size dependencies are not stated ex- 
plicitely. In what follows, all energies are measured in units of the Ising coupling constant 
J, all Temperatures in units of J/ks {kB'- Boltzmann's constant). 

B. Some remarks on Monte-Carlo simulations 

For a Monte Carlo simulation, a Markov process is set up on configuration space Tj^d 
with a certain problem adapted stationary distribution w{S), which is assumed^ to depend 
only on the interaction energy E and the magnetization M of the microstate S, i.e. w{S) = 
w{E{S),M{S)). From the Markov chain {S}j^ of length J\f (which, at least in the limit of 



^It is straightforward to extend the formalism introduced in the following sections to the case of a 
more general stationary distribution. The restricting assumption w{S) = w{E{S),M{S)) is made 
only for the sake of notational simplicity. 



infinitely long samples, is distributed according to w) the simulation average of an arbitrary 
function / : rj;^d ^ R on configuration space is obtained: 

(/(^)U.({'5W) ■■=^11 /(^) ""^ E fis)HS). (2) 

Of course, the simulation average depends on the stationary distribution and, unless the 
length of the Markov chain reaches infinity, on the particular sample {S}j^. 

C. Standard histogram technique 

The histogram H^{E, M; {S}j^), which is proportional to the number of microstates 
of the sample with interaction energy E and magnetization M, is given by the simulation 
average of the observable Se{s),eSm{s),m'- 

T? z_^ Se(s),eSm{s),m 



5G{5U 



N- 



E '^S(5),E'^M(5),mw(E(S'),M(S')) 

= n{E,M)w{E,M) . (3) 

Since the underlying stationary distribution w is known (at least beside an irrelevant factor), 
the density of states is obtained as 

Q{E,M) := 2_^ ^e{s),eSm{s),m 

= hm HUE,M-{SU)/w{E,M), (4) 

or — more realistically for a computer simulation — at least an estimator for ^2 is obtained 
by omitting the limiting procedure lim^v'^oo- 



D. Transition observable method 

In this section it is shown that the density of states can be obtained from simulation 
averages of certain transition observables defined below, which have the advantageous feature 
that they enable the estimation of the density of states in a much more efficient way than 
the standard histogram method does. 

As a preliminary step, let us define the microcanonical average of any system observable 

fisy. 

{m)iE,M):= hm (/(^)^^(^)^^W).m) ({^M 

^^^ ^^^ ^ ^— HUE,M;{SU) 

Yl ^e{s),eSm(s),m fi^) 

Q{E,M) ■ ^ ' 



Let ^ be a set of operators acting on configuration space F^d 

AC {A: AS eVL-i y S e 1^4 . (6) 

The transition observable Nj^{S) is defined as the number of operators A E A acting on 
the particular microstate S, which result in microstates S with interaction energy E{S) = 
E{S) + iAE and magnetization M{S) = M{S) +jAM: 

^A (^) -^ Z^ ^E(S),E{S)+i-AEKl{S),M{S)+j-AM ^"AS,S '■> %■,] El- ■ (7) 

(See Fig. |l] for an illustration of the thus defined observables.) Then, for any set of operators 
A which satisfies 

7^ 2_^ Se{S),eSm{S),M ^A^^) = Z_^ SE{S),E+i-AESM{S),M+J-^M ^a'' \S) y (§) 

the density of states Q{E, M) can be calculated from the microcanonical averages (^ of the 
thus defined transition observables to yield 

(tV^* "•''(5)) {E+iAE,M+jAM) 

Q{E,M) = Q{E+iAE,M+jAM) . (9) 

{N'j^'(S))iE,M) 

Remarks: 

1. Microreversibility as explained in App. ^ is a sufficient condition for the equality in 
(^. Apart from this condition which is implemented easily, the set A of operators can 
be chosen arbitrarily. 

2. In the density of states, a multiplicative constant is physically irrelevant! For that 
reason, Q can be chosen arbitrarily for one particular value of {E,M). Then, the 
density of states of the remaining [E, M) values can be calculated from Eqn. (P). 

3. The efficiency of the transition observable method depends crucially on the particular 
choice of A. 

4. The transition observable method is neither restricted to the investigation of Ising 
systems (with bare next neighbour interaction, cf. |T^) nor to the investigation of dis- 
crete systems (cf. [|^). Example: consider a discrete spin system with a Hamiltonian 
consisting of two interaction terms 

n{S) = E,{S) + E2{S), (10) 

which depends on certain coupling constants, say, Ji and J2 (e.g. ferromagnetic cou- 
pling to next neighbours and antiferromagnetic coupling to next-nearest neighbours). 
The knowledge of the density of states as a function of Ei and E2, i.e. 

n{Ei,E2) := ^ 5E^{S),EiSE2iS),E2 ■ (11) 



enables the determination of the thermostatic properties of the system for all possible 
values of the ratio of the coupling constants by applying certain " skew- summing" 
techniques (cf. [|13)- In complete analogy to the above, a set of transition observables 
can be defined, which facilitates the determination of the thus defined density of states 

5. For a matrix T defined as 



1 / e'-e m'-m \ 

m,., ,.,w. ,n := ^^(a^^^"^ ' ^" iS)jiE,M) , (12) 



liE',M'),iE,M) ■- 1^1 y 



where |^| is the cardinality of the set A, it is easy to show that 

a) T is a stochastic matrix, i.e. 

miE'M'UEM) ^ ^ (^'' ^'), (E, M) (13) 

and 

/ . V^\{E',M'),{E,M) ^ ^ ^^^^ 

(E',M') 

b) The density of states is the stationary state of T: 

E [T](E',M'),(i.,M) ^(E, M) = n{E', M') . (15) 

(E,M) 

Furthermore, if T is regular, i.e. if for all £", M' and E^ M 

3n e N : [T'^\e'M'Uem) > Vm > n , (16) 

the stationary state of T is unique. Nevertheless, even if this condition is not fulfilled, 
the density of states can be computed piecewise on certain subsets £,M. of the total 
set of possible energy and magnetization values of the system. The thus produced 
"fragments" of the density of states are then connected to each other via (a priori 
unknown) multiplicative constants. 
Note here that the stochastic matrix T defined above is closely related to the so-called 



Transition Matrix Monte Carlo method introduced by Swendsen and Li [|T5|. 



E. Reduction to Oliveira's observable: the reduced transition observable method 



The results of Sec. |11D| can be simplified to those presented by Oliveira in reference [Q, 



where no information on the magnetization of the system is regarded. Formally, this can be 
achieved by a summation over the magnetization M (or the index j, respectively) in some of 
the expressions of the preceding section. Then, however, only a determination of the reduced 
density of states 
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f2(E) :=^fi(E,M) (17) 

M 

is feasible, which does not entail the entire thermodynamic information of the system (in 

the sense that Q{E, M) enables the calculation of thermal and magnetic equations of state 

in various ensembles whereas Q{E) just allows for the estimation of the thermal equation of 

state). 

The reduced microcanonical average of any system observable f{S) over the energy-shell 

E{S) = Eis defined: 

^^^'^^ ^^^ = ^'^^ EHUE,M;{SU) = m • (^«) 

M 

We further define the reduced transition observable 

= /_^ "E(S),E{S)+i-AE /_^ ^AS,S 5 «, J G Z . (19) 

5er,d A&A 



L" 



Then, for any set of operators A which satisfies 

^ Y. ^EiS),EN\iS) = Yl ^EiS),E+iAEN^\S) , (20) 

the reduced density of states ( p^Tf ) can be calculated from the reduced microcanonical average 



of the reduced transition observable (|T^) 



{N'^{S)){E) 

Q{E+iAE) = Q(E) . (21) 

{N;^\S)){E+iAE) 

Again, microreversibility (cf. appendix 0) is sufficient to ensure the equality in (|20|) . 



III. COMPARISON OF THE EFFICIENCY OF THE STANDARD HISTOGRAM 
AND THE TRANSITION OBSERVABLE TECHNIQUE 

To demonstrate the advantages of the transition observable method, numerical results 
obtained from either the standard histogram method or the transition observable method 
are compared. The philosophy of the comparison is to use some simulation technique to 
generate one sample of microstates which then is evaluated according to both methods. 
The simulations were performed for a d = 2, L = 32 and a d = 3, L = 10 Ising system. For 
the sake of completeness, the details of the computer simulations are given in App. y. 
The set A of lattice operators was chosen to consist of L'^ operators, which are labelled by 
the subscript i and are defined by their action on a particular microstate S: 



Ai : S = ai,a2,...,(Ti,...,aLd 



S = ai,a2, .-., —ai, ..., aid 



(22) 



i.e. the operator Ai flips only the i-th spin of the Ising lattice. Obviously, since AiAiS = S, 
the thus defined set of operators meets the condition ( |A2| ) and therefore is microreversible. 
Note that the determination of the simulation average of the transition observable by use 
of this particular set of lattice operators can be done very fast. In fact, the time needed for 
applying the L"^ operators of the set .4. to a particular microstate is much shorter than the 
time needed to perform a lattice sweep! 

Simulation averages of A''^-' were recorded only for values ofi G { — 1,0,1} and j G 

{-1,1}- 

To emphasize the difference between the two methods, we compare "discrete derivatives" 

(i.e. ratios of differences) of the logarithm of the density of states, namely: 



(i) For the case oi d = 2 Ising model: 



As(^lnfi(i?)) : = 



1 

2AE 



-^— In 



\n^{E+AE) — IndiE-AE) 

{n\{S)){E-AE){n\{S)){E) 



{N;^\S))(E+AE){N^\S))iE) 



(23) 

(24) 



(ii) For the case of d = 3 Ising model: 



Am ( \nQ.(E,M) 



1 



2AM 



1 



"2 AM 



1 



■2AM 



1 



■2AM 



In 



In 



In 



In r2(£;,M+AM) — In r2(_B,M-AM) 



{n^^{S)){E,M-AM) (^N^^''^{S)){E+AE,M) 



(iVJj'"^(5))(£,M+AM) 



'{S)){E+AE,M) 



{N-^-^{S)){E,M-AM) (^N^^{S)){E~AE,M) 
{n^'^'''^{S)){E,M+AM) {n^''^{S)){E~AE,M) 



{n°/{S)){E,M-AM) {N°/{S))iE,M) 
[N^-^{S))iE,M+AM) {n^~'^(S))(E,M) 



(25) 



(26) 



(27) 



(28) 



Note that both "discrete derivatives" (|23| ) and ( pSj) are closely related to microcanonical 
equations of state (see ||^11| and appendix^ for more details). 



A. Example 1: the d = 2, L = 32 Ising lattice 



In Fig. ^, the differences of the logarithm of the reduced density of states as emerging 
from the transition observable method (cf. (p^ ) and the conventional histogram method 
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(cf. (^) are shown together with the exact result |0. By use of a sample of 8 ■ 10^ mi 



crostates, the transition observable method yields a result which, on the scale of the figure, 
can hardly be distinguished from the exact result, whereas the data obtained from the his- 
togram method scatter strongly around the latter. 

In Fig. 1^, the results of the transition observable method (using one sample of n = 8 ■ 10^ 
microstates) are compared to the results of the histogram method for several sample lengths 
(n, 5 • n, 10 ■ n and 15 ■ n) by plotting the deviation of the simulation data from the exact 
result. Even if the simulation time is chosen 15 times longer in the histogram method, the 
transition observable method still yields more accurate results. 

Calculating the mean square deviation of the simulation data from the exact result as a func- 
tion of simulation timeQ (Fig. ^, we notice that the accuracy of both methods is improved 
according to a power law (the corresponding exponents seem to be the same (~ —1) in both 
methods). But, at any given time, the transition observable method beats the standard 
histogram method in accuracy by a factor of roughly 40. 



B. Example 2: the d = 3, L = 10 Ising lattice 

From the Monte Carlo samples, we computed the differences of the logarithm of the 
density of states in direction of the magnetization according to Eq. (^ in the case of the 
histogram method and according to Eqs. (p6|)-(p8l) in the case of the transition observable 
method (in fact, we computed the mean value of the three possibilities (|26D-(|28D of deter- 
mining Am (Infi)). 

In Figs. I^a) and b), A^/ (In Vt) is shown for E /IQ^ = —.924. In Fig. ^), a sample of length 
10-10^ microstates is used for the evaluation of A^ (In fi) according to both methods whereas 
in Fig. ^), the histogram method with a sample length of 50-10^ microstates is compared to 
the transition observable method with sample length 10-10^ again. For a better visualization 
of the difference between the two methods, an odd polynomial (ffu^M) = a-M+b-M^+c-M^) 
was fitted to the transition observable data. Subtraction of the data of Figs. ^) and b) from 
this polynomial yields the plots shown in Figs. §c) and d). In the figures, the transition 
observable (histogram) data are represented by the solid lines (points). The plots of the 
differences show the consistency of both methods, that is: both data sets scatter "randomly" 
around the fit function. The data emerging from the transition observable method, however, 
are much more accurate than the data emerging from the standard histogram method even 
if much longer samples are used in the latter. 

For a quantitative comparison between the two methods, the mean square deviations 
of the simulation results, with respect to a fitQ to the best data obtained by the transition 



^The comparison of the data was done within a certain "energy window" which was chosen around 
the centre of the histogram, i.e. the tails of the histogram have been discarded. Since the same 
sample is used in both evaluation techniques, the result of the comparison does not depend on the 
width of the "energy window" chosen for the evaluation of the x'^-deviations. 

^We have performed a weighted x^ ht of an odd polynomial ffn^M) = a- M + h ■ M^ + c • M^ to 



observable method, is shown as a function of the simulation time in Fig. |^. The accuracy 
is improved according to a power law with exponent ?« — 1 in both methods. But: at any 
given time, the transition observable method yields results which are more accurate than the 
results emerging from the histogram method by a factor ^ 1/250 in the sense of the mean 
square deviation. In order to produce results of similar quality, the simulation time in the 
standard histogram method has to be ~ 250 times longer than in the transition observable 
method! 



C. General remarks on section III 



1. The two examples discussed in the preceding sections show that a simulation can be 
accelerated dramatically by use of the transition observable method. Here, it is not 
the algorithm to speed up the simulation but it is the observable measured during 
the simulation! The reason for this striking difference is indeed very simple: while 
every microstate, which is decided to be part of the sample, just yields one entry 
in a list in the histogram method, it might yield many transitions to neighbouring 
states (neighbouring with respect to the interaction energy and/or magnetization) and, 
hence, the statistics of the transition observable method can be expected to be much 
better than the statistics of the conventional histogram method. In fact, since the Ising 
systems under consideration just allows for 5 (7) different interaction energy changes 
and only two magnetization changes^ under single spin flip operations [in d = 2 (3)), 
the set of operators A chosen above can be expected to shorten the computational 
effort by a remarkable factor, roughly proportional to the square-root of the inverse 
volume L~'^/^ of the system! 

2. The change of interaction energy under a single spin flip operation depends on the 
configuration of the spins in the very neighbourhood of the particular spin to be 
flipped. The typical configurations of neighbouring spins vary with the interaction 
energy of the whole system. For that reason, the factor of proportionality of the 
efficiency gain in the sense of the x^-comparison introduced above can be expected to 
depend on the mean interaction energy of the histogram, which itself depends on the 
simulation parameters (i.e. the stationary distribution). 

3. The particular way of generating the sample of microstates is not important in the 
context of the comparison of the two methods introduced in Sec. ||. 



the data obtained from a 50 • 10^ sample by applying the transition observable method. The errors 
needed for the weighted fit have been produced by a jack-knife blocking procedure using 25 data 
sets of length 2 • 10^ sweeps. 

^AE/J = (±12), ±8, ±4, and AM = ±2 in d = (3), 2. 
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IV. CONCLUSION 

A Monte Carlo simulation consists of two steps. The first step is the generation of a sam- 
ple or spot check of microstates. The second step is the investigation of these microstates. 
Conventionally, if the aim is to speed-up the simulation, the first step is modified while 
the second remains unchanged. We have shown that a more extensive exploitation of the 
microstates of the sample, i.e. taking simulation averages of the transition observable in- 
stead of just cumulating a standard histogram, can effectively speed-up the simulation by a 
tremendous amount! 

In an extremely straightforward implementation of the transition observable method, 
we reach a speed-up which can be expected to be proportional to the square root of the 
volume L'^ of the system under consideration. Such a speed-up seems unattainable by an 
improvement of the algorithm of the simulation, i.e. by modifying the first step of the 
simulation. 

Even though the transition observable method seems to be built for discrete spin systems, 
one of us (J.D. Munoz, cf. |jl3|) has already shown that the method can be transferred to 
continuous spin systems. 

An extension of this method to "non-spin" systems like polymers might be a topic of 
future investigations. 
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APPENDIX A: MICROREVERSIBILITY 

Let ^ be a set of operators acting on configuration space V^d 

A(^{A: ASeVLdM S eVLd} (Al) 

such that for all A E A, there exists a unique inverse operator B = A~^ G A, i.e. 

\fAeA 3! Be A: BAS = S . (A2) 

Then, A is said to show microreversibility. 

From the microreversibility of A, it follows immediately that the number of operators A e A 
which transforms 5* into 5* equals the number of operators A e A which transform 5* back 
into S, i.e. 

Yl ^AS,S = Yl ^s,As ■ (A3) 

AgA AeA 

Using the definition of the transition observable (0) , this can be shown to be equivalent to 
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/ ^ Se(S),eSm{S),M ^a{^) = Z_^ SE{S),E+i-AESMiS),M+j-AMN_/' \S) . (A4) 

That is, the number of operations which transform microstates with interaction energy E 
and magnetization M into microstates with E + iAE and M + jAM by use of operators 
A & A is identical to the number of operations which transform "backwards", i.e. from 
states with E + iAE and M + jAM to those with interaction energy E and magnetization 

M. 



APPENDIX B: MICROCANONICAL EQUATIONS OF STATE 



As mentioned in Sees. |T|, the differences of the logarithm of the density of states 
Am (In Q) are related to the microcanonical magnetic equation of state. Indeed, equation 
(^) is the microcanonical magnetic equation of state in a discrete notation (appropriate for 
the description of finite Ising systems), which converges in the thermodynamic limit L — > cxo 
towards the magnetic equation of state of the infinite system 

-^{e,m) = ^ hm L-'^lnfi(E, M, Z.-^) , (Bl) 

where e := L~'^E and m := L~'^M are intensive quantities. The difference of the logarithm 
of the reduced density of states, as defined in Eq. (^31), converges towards (3{e, h/T = 0) of 



the infinite system for zero external field and can serve to compute zero field properties of 
the system^. 

Note that it is unnecessary and a rather roundabout way to convert the thus obtained 
data into the commonly used canonical quantities. For details on the investigation of phase 
transitions in a microcanonical approach and a microcanonical finite-size scaling theory see 
references |Q , ||T^ and |Tl[ . 



APPENDIX C: DETAILS OF THE MONTE CARLO SIMULATION 

1. Simulation of the d = 2, L = 32 Ising lattice 

A Monte Carlo simulation of a 32^ Ising system with periodic boundary conditions was 
performed. The stationary distribution of the underlying Markov-process was chosen to be 
proportional to the Boltzmann weight w{E{S),M{S)) oc exp {—TC{S)/T} with simulation 
parameters h = and T = 2.269 (see Sec. P A| for the definition of the Ising Hamiltonian). 
We have implemented a sequential lattice update with a "Metropohs-type" transition rate 
T(5'— >5") = m.m{l, w{S')/w{S)} and we have sampled every L^ configuration only. After 
"equilibration" (6.4 ■ 10^ lattice sweeps have been discarded), several successive samples of 
8-10^ microstates were taken. 



^/3(e, h/T) is the derivative with respect to e of the Legendre transform of 
linii^oo ^ '^lrir2(e, ?n,, L"-*^) with respect to m. 
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2. Simulation of the d = 3, L = 10 Ising lattice 

A Monte Carlo simulation of a 10^-Ising system with periodic boundary conditions 
was performed. The stationary distribution was chosen to be w{E{S),M{S)) oc {{Eq — 
E(S)) / Nq}'^^°~'^^^'^ , i.e. independent of M{S) again. The parameters were chosen to be 
Eq = 1586 and A^o = 1000 (for a detailed discussion and interpretation of this stationary 
distribution, see [|l9l). The way of updating the lattice configurations is the same as for the 
simulation of the 32^-Ising system (cf. App. |C 1| ) . After "equilibration" (2-10^ lattice sweeps 
have been discarded), several successive samples of 2 • 10^ microstates were taken. 
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Figure 2: 
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Figure 3: 






<1 



0.004 - 

0.002 - 

- 

-0.002 - 

-0.004 - 

-1.8 



8 ■ 10^-sweeps 
8 ■ 10^-sweeps 

T r 



T 




-1.6 -1.4 -1.2 









8 ■ 10^-sweeps 
5 ■ 8 • 10^-sweeps 

T 




-1.6 -1.4 -1.2 
E/32^ 









8 ■ 10^-sweeps 
10 • 8 ■ 10^-sweeps 

T 




-1.6 -1.4 -1.2 

E/322 









8 • 10^-sweeps 
15 ■ 8 ■ 10^-sweeps 

T 




-1.6 -1.4 -1.2 

E/32^ 



Figure 4: 



-7p 

-8 - 

-9 - 

-10 - 

-11- 

^J -12 - 

^ -13 - 

-14 - 

-15 - 

-16 - 

-17 







X^-deviation from the exact result 



_L 



1 ~ — I X~~ 

mstogram method 

transition observable method 




X " X X XX, 




_L 



_L 



_L 



12 3 4 

ln{simulation time) 



19 



Figure 5: 
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Figure 6: 
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Captions: 



FIG. 1. a) Visualization of the transition variables N^J^ for the case i G {—1,0,1} and 
j € {—1,1}. Given a particular microstate S with interaction energy E{S) = E and magneti- 
zation M{S) = M, Nj^\S) gives the number of possibilities to reach any state S with energy 
E{S) = E{S) + i ■ AE and magnetization M{S) = M{S) + j ■ AM by applying the set of lattice 
operators A to the microstate S. The microcanonical average of the transition observable A^V is 
proportional to the total number of possibilities for the event that, given any state S with energy 
E{S) = E and magnetization M{S) = M, any other state S with energy E{S) = E{S) +i- AE and 
magnetization M{S) = M{S) +j ■ AM is reached under the action of A. Fig. ||b) shows the "tran- 
sition paths" corresponding to the various microcanonical expectation values contributing to the 
differences of the logarithm of the density of states Am (InO) at the point {E,M), as introduced 



in Sec. III. 



FIG. 2. Intensive energy as a function of 



-1 



of a 32^ square Ising system. The 



Aij(^lnOj 

solid line is the exact result, the dashed line corresponds to the transition observable method and 
the points to the conventional histogram method. The same sample of 8 ■ 10^ microstates was 



used to perform both evaluations. The energy is plotted against A^; (in $7) because of its 
correspondence to a thermal equation of state E{T); cf. App. ^. The strong fluctuations in both 
the high energy and the low energy region of the figure are due to poor statistics in the tails of 
the histograms. In the central region, the difference between the dashed and the full line is smaller 
than the line thickness! 



FIG. 3. In order to point out the differences between the transition observable method and the 
conventional histogram method and in order to show that the used estimators are indeed unbiased, 
the data emerging from the Monte Carlo simulation are subtracted from the exact result. The 
transition observable data (histogram data) are represented by solid lines (points). 



FIG. 4. In order to judge the quality of the two methods, the x^-deviation of the simulation 
results from the exact result is shown as a function of the simulation time (in units of 8 • 10^ lattice 
sweeps). 



FIG. 5. a) and b): differences of the logarithm of the density of states Ajv/(lnO) for 
-E/10^ = —.924. In Fig. a), a sample of length 10 • 10^ microstates is used for the evaluation 
of Am (Inil) according to both methods whereas in Fig. b), the histogram method with a sample 
length of 50 • 10^ microstates is compared to the transition observable method with sample length 
10 • 10^ again. For a better demonstration of the difference of the two methods, the same data are 
subtracted from a fit function in Fig. c) and d) (see text for the details of the fit). In all figures, 
the transition observable data (histogram data) are represented by the solid lines (points). 
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FIG. 6. In order to compare the quality of the results emerging from the histogram method to 
those emerging from the transition observable method, the mean square deviations of the simulation 
results with respect to a fit to the best data obtained by the transition observable method is shown 
as a function of the simulation time (in units of 2 • 10^ lattice sweeps). 
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